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Abstract 


A theory  of  a-posteriori  estdjnates  for  the  finite  element  method  was 
developed  earli'ar  by  the  anthors.  Based  on  this  theory,  for  a two-po^Jit 
boundary  value  problem  the  existence  of  a unique  optimal  mesh  distribution 
is  proved  and  its  properties  analyzed.  This  mesh  is  characterized  in 
terms  of  certain,  easily  computable  local  error  indicators  which  in  turn 
allow  for  a simple  adaptive  construction  of  the  mesh  and  also  permit  the 
computation  of  a very  effective  a-posteriori  error  bound.  While  the  error 
estimates  are  asymptotic  in  nature,  numerical  experiments  show  the  results 
to  be  excellent  already  for  10%  accuracy.  The  approaches  are  not  re- 
stricted to  the  model  problem  considered  here  only  for  clarity;  in  fact, 
they  allow  for  rath  r straightforward  extensions  to  more  general  problems 
in  one  dimension  as  v?ell  as  to  higher  order  elements. 


1.  Introduction 


For  the  numerical  solution  of  boundary  value  problems  by  finite-element 
techniques t the  construction  of  optimal,  or  near-optimal  meshes  Is  of  con- 
siderable practical  importance.  The  same  can  be  said  when  finite-difference 
or  collocation  methods  are  used.  Many  articles  in  the  literature  deal  with 
questions  that  bear  a relation  to  this  problem,  yet,  as  observed  in  Cl5], 
even  for  two-point  boundaz*y  value  problems  relatively  few  address  it  directly. 
We  shall  not  atten^Jt  to  survey  this  literature. 

There  are  various  ^U1alyses  of  the  approximation  error  of  a given  func- 
tioj  by  piecewise  polynomials  with  a fixed  number  of  pieces  of  fixed  order 
(see,  e.g.,  [7],  [8],  CIO],  [12],  [21]  and  the  references  cited  there).  In 
principle,  such  studies  may  relate  to  the  finite  element  method  since  that 
method  leads  to  optimal  igiproximations  under  the  energy  norm.  The  mentioned 
error  estimates  involve  higher  derivatives  of  the  given  function.  With 
these  results  as  a basis,  a nundser  of  authors  developed  methods  for  the 
construction  of  optimal  meshes  for  ooUocation  and  finite-difference  methods 
([11],  [13],  [14],  [17],  [22],  [23]).  For  this  the  needed  information  about 
the  derivatives  of  the  solution  is  obteuLned  from  the  approxinate  solution, 
for  instance,  by  means  of  difference  fomilas.  This  procedure  can  be  theo- 
retically justified  in  the  case  of  r^ular  meshes  (see,  e.g.,  [20]).  How- 
ever, When  there  am  abnipt  changes  in  the  mesh— as  they  arise  with  refine- 
ment techniques— then  "internal  boundary  l^^rers"  appear  in  the  error  function 
(see»  e.g.,  [2],  [6]),  and  hence  the  difference  formulas  cease  to  approximate 
well  the  desired  derivatives. 
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Various  results  on  optimizing  finite  element  meshes  have  appeared  in 
the  engineering  literature.  Without  entering  into  any  detedlsy  we  mention, 
for  Instance,  the  articles  C9],  [18],  [19],  [25],  [26],  and  [27]. 

In  recent  year,  for  initial  value  problems  for  ordinary  differential 
equations  very  effective  procedures  have  been  designed  and  arudyzed  for 
adapting  the  stepsizes  and  the  order  of  the  numerical  methods  (see,  e.g., 
the  survey  [IS]).  The  principal  tool  is  the  availability  of  error 
analysis  with  a local,  a-posteriori  character.  These  estimates  are 
asymptotic  in  nature t yet  practical  experience  has  proved  their  relia- 
bility for  reasonable  tolerances. 

In  this  paper  we  use  a new  approach  to  the  construction  of  optimal 
finite  element  meshes.  It  is  based  on  a theory  of  a-posterlori  estimates 
for  the  finite  element  method  developed  in  [3],  [4]  (see  also  [$]).  As 
in  the  case  of  the  initial  v^due  problems,  the  estimates  are  asyngxtotic 
in  character.  More  specifically,  higher-order  terms  in  the  maximal  mesh- 
size  H are  neglected;  that  is,  asyndetic  expressions  of  the  form 
1 ■*’  o(l)  as  0 are  considered  to  be  approximately  equal  to  one.  At  the 
same  time,  edl  constant  factors  of  these  (l-i'o(l))-’terms  can  be  evaluated 
computationally. 

For  clarity  of  presentation,  we  restrict  the  discussion  to  a sijiq^le 
two-point  boundary  value  problem  involving  a linear,  self-adjoint  operator 
of  second  order.  Moreover,  for  sin^licity,  we  erploy  only  piecewise  linear 
elements.  The  approaches  allow  for  rather  stralghtforMBrd  extensions  to  a 
variety  of  more  general  problems  in  one  dimension,  and  there  are  no  essen- 
tial limitations  to  the  use  of  higher-order  elements.  In  fact,  analogous 
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techniques  even  permit  consideration  of  elements  of  different  order  in 
different  parts  of  the  mesh* 

Continuous  mesh  distribution  functions  are  used  to  prove  the  existence 
of  a unique » optimal  mesh  distribution  and  to  analyze  its  properties.  In 
particular,  it  is  shown  that  the  value  of  the  optimal  error  is  rather 
stable  under  perturbations  of  the  optimal  mesh.  Hence  it  is  indeed  un- 
necessary to  compute  this  mesh  with  excessive  accuracy.  The  optimal  mesh 
is  characterized  in  terns  of  certain  easily  computable  local  error  indica- 
tors. This  allows  for  a single  adaptive  method  to  oonstruct  that  mesh 
(see  [4])  and,  at  the  same  time,  to  confute  very  effective  a-posteriori 
error  bounds.  Although,  as  mentioned,  the  error  estimates  have  only 
asyndetic  character,  numerical  experience  shows  that,  as  in  the  case  of 
initial  value  problems,  the  results  are  excellent  already  for  aoouraoies 
of  the  order  of  ten  percent* 
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2.  Notation 

Let  I z I(a,6)j  a < 3,  be  the  open  interval  (x  € R^;  a < x < 3)  and 
I = i(a,3)  its  closure.  As  usual,  H^(I)  denotes  the  space  of  square- 
integrable  functions  on  I and  C^(I)  c H^(I)  the  subspace  of  continuous 
functions  on  I.  The  norms  on  H^(I)  and  C^(I)  will  be  written  as 
j||*||q  and  respectively. 

Define  E(I)  as  the  space  of  real,  infinitely  differentiable  functions 
on  I for  which  all  derivatives  have  continuous  extensions  on  1.  Mcreover, 
let  V(I)  c E(I)  be  the  subspace  of  all  functions  with  compact  support  in  I. 
For  any  integer  k > 1,  the  spaces  H^(I)  and  C^(I)  are  the  conpletions  of 
E(I)  under  the  norms 


(2.1)  I 

and 


(2.2) 


> 


respectively.  Analogously,  the  completions  of  P(I)  under  these  twc  norms 
are  the  spaces  H^(I)  and  c^d). 

Let  a,b  € C^(I)  be  given  such  that  a(x)  ^ a > 0,  b(x)  >0,  V x ^ I. 
Then  E(I)  and  E^^d)  shall  be  the  spaces  H^(I)  and  H^d),  respectively. 
with  their  norm  replaced  by 
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(2.3)  j||u|||  = [ / (au*^+bu^)dx],  (u'  = . 

If  b = 0 on  I,  then  (2.3)  is  only  a senunom  on  E(I).  On  the  other  hand, 
on  Eq(I),  (2.3)  is  alvjays  a nom  which,  noreover,  is  equivalent  to 
Obviously,  Eq(I)  is  a Hilbert  space  and  for  b ^ 0 the  same  is  true  for 
Ed).  We  denote  the  inner  product  in  either  space  by  j(.,.)^.  For  b = 0 
on  I,  E(I)  is  a Hilbert  space  modulo  the  constant  functions. 

On  I = i(a,B)  we  consider  partitions 

(2.4)  A(' ) : a = Xq  < Xt  < . . . < < x^  = B,  m = m(A)  5 1, 

and  introduce  the  notations 

Ij(A)  = I(x^_3^,x4)' 

► 3 “ ,m 

(2.5)  h.(A)  = x^  - x^  , 

3 J j-J-  j 

h(A)  = max  h.(A),  h(A)  = min  h.(A) 

j ••1 5 ♦ • • jin  j ••X  5 • • • jin 

All  partitions  A which  for  fixed  X > 0,  x > 1 satisfy 

J (2.6)  h(A)  > XMA)*^ 

f are  said  to  be  (X,k) -regular. 

For  given  A = A(I),  we  denote  by  S(I,^)  c H^(I)  and  Sq(I.A)  c hJ(I)  the 
sui^Si-aces  of  all  functions  for  which  the  restriction  to  any  3 = l,...,m, 

is  linear.  Analogously  P^(I,A)  c H^(i)  and  P^(I,A)  c Hg(I),  k > 0,  shall  con- 
sist of  the  functions  for  which  the  restrictions  to  Tj(A),  3 = l,...,m,  belong 
to  C^(Ij). 
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For  later  use  we  lote  tte  following  well-known  leniaa  (see,  e.g. , [2*1]} : 
Lenca  2.1;  For  given  I = l(a,0),  a < 8,  and  A = A(I)  there  exists  a posi.- 
tive  constant  K such  that 


(2-7)  inf  < KWA)^)^-,  V u € H^(I)  0 1^(1) 

w€SQfl[,A)  I E I 2 *0 
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3.  A Boundary  Value  ProbLec 
3.1  Basic  FonHulation 

As  mentioned  in  the  introduction,  we  restrict  the  discussion  to  a 
ginipio  nrviol  probleQ.  FoT  ease  of  notation,  the  unit  interval  I = 1(0,1) 
is  used  from  now  on  throu^iout  the  remainder  of  Ihe  ^aper.  On  I we 
consider  the  equation 

(3.1)  Liu]  = - ^ a(x)  ^ + b(x)u  = f(x),  X € 1(0,1)  , 
together  with  the  boundary  conditions 

(3.2)  u(0)  = u(l)  = 0 . 

We  assume  that  a € C^(I),  b,f  € C^(I),  and,  as  before,  that  a(x)  > a > 0, 
b(x)  ) 0,  V X € I. 

The  weak  solution  of  the  problem  is  the  unique  Uq  € Eq(I)  with 

(3.3)  ^ ^ ^ 

where 

(3.4)  F^(v)  = / fvdx  . 

I 

Note  that  under  our  differentiability  assunptions  about  a,b,f  the  solution 

3 

Uq  of  (3.3)  belongs  to  C (I)  and  also  satisfies  (3.1/2). 

With  the  partition  A and  the  space  Sq(I,A)  c H^(I)  specified  as  in 
Section  2,  we  consider  the  finite  element  solution  u^  € Sq(I,A)  defined  by 


(3.5) 


i(u^,v)e  = F^(v),  V V € Sp(I,A)  - 


Since  u^,  € H^(I)  n 1^(1),  it  then  follows  ftorn  Lenina  2.1  that 

(3.6)  ^ ^ 5'^^A)ji|uQ||2  - 
3.2  A-Ppsteriori  Error  Analysis 

ffe  ocHisider  the  residual  r = L(u.)-f  on  the  intervals  I.,  that  is, 

a J 

(3.7)  r-(x)  = (L(u.)-'f)(x),  V x € I.,  j = l,...,m  . 

3 A 3 

Let  z.  X EqCI^)  be  the  solutions  of 

(3.8)  T (2-,v)p  = F (v),  V V € E^(I.),  3 = l,...,ni 

ij  J £.  Tj  3 

and  set 

(3.9)  Z(l)^  = Z T • 

3=1  ^3  ’ ^ 

The  following  result  was  proved  in  [ 3 ] : 

Theorem  3.1;  The  error  e = u^  - Ug  satisfies 

(3.10)  jllellE  = Z(5)^(l+0(h(A)))  as  h(A)  0 , 

where  the  (constant  in  the  bound  of  the  O-tem  depends  on  a and  b 
but  not  on  f and  A. 

We  analyze  the  quantity  Z(A)  of  (3.9)  in  sane  more  detail.  Because 
u^  is  linear  on  any  I,,  we  nay  write 

rj  = L(u^)-f  = -a'u^  t bu^  - f r , 


(3.11) 
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Hhe9% 

(3.12) 

(3.13) 


p-(x)  = a(x)Up(x),  T.(x)  = -a*(x)e'(x)  + b(x)e(x),  V x € I-  . 
J J 3 

€ Eg(Ij)  be  such  that 

(a)  j (<p.,v)  = (v),  V V € E.(I.) 

Ij  ] t p.  0 3 

(b)  i/tj»v)j.  = (V),  V V € Epdj)  , 


and  therefcjre  z.  = (pj  'kj* 

The  smallest  eigenvalue  of  the  differential  operator  L on  with 

2 2 

zero  boundary  conditions  is  bounded  below  by  the  smallest  eigenvalue  an  /h^ 

2 2 

of  the  operator  /dx  on  I^.  Hence  it  follows  from  (3.13b)  that 

iJl'kjilo  ^ Chj 

1 • 1 


and  therefore 


(3.14) 


(>|rj)  ^ 


2 

0 * 


Here  as  in  subsequent  estimates  C denotes  a generic  constant  which 
has  different  values  in  each  instance  but  is  independent  of  the  other  essen- 
tial variables  in  the  same  expression.  Now  note  that 


(3.15)  ^ = I (-a'e'+be)^dx  ( C J [(a'e’ )^+(be)^]dx  ^ C , ||e||p  , 

^ “ Ij  Ij  ^ 

which  together  with  (3.14)  gives 


.If. 


(3.16)  ^ 

Ue  iMindBioe  (g^iasislities 


aad 


Hi}  = 


" m 1/2 


(3.17) 


Rcaa  Thecsran  3.1  ue  obtain  — fdth  soae  fo|  5 1 — 


jlei^  = AA)(l+0(h))  = 


m 


I (^j^j  *#j-»t-)£(l+0(h)) 


= [Q(A)^*2&Q(A)jieS£R(A)+jie|^(A)^](l+<Xh)) 


= (Q(A)+aj||ej|j^(A))^(l40(h))  + (l-o^)j||e!||R(A)^(l-K)(h)) 

But  by  (3.16)  v#e  have  R(A)  = 0(h)  and  thus 

jilell^  = (Q(A)+trj||e|ij^(A))^(l+0(h))  , 
vdiich  in  tum  irplies  that 
(3.20)  jllell^  = Q(A)(l+0(h))  . 

Therefore,  in  view  of  (3.10),  we  have  proved  the  following  result: 
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ThBPgwBi  3-2:  Qrs  lisi)  aodl  <^ii}  Sie  <Sefisw&  by  (3.S)  acad  (3.15j>,  srspeo^ 
niaidly- 

(3.21)  QW^  = Z(A)^(1-K)m)>,  as  Ma)  - ®, 

«itez«  the  oonsttanit  iss  iSie  boaoid  of  the  O-ttm  oippPBwte  on  a ara!  b but 
not  on  f and  A. 

iiiile  ^pmadnations  of  Z(A)  can  be  oonputed  and  Q(a)  is  not  re^lily 
accessible,  the  quantit:,  Q is  better  suited  than  Z for  our  theoretical 
studies  of  optical  partitions  A. 
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4,  Optiaal  Partitions  - Case  I 

In  this  section  we  restrict  ourselves  to  the  case  tAien  ^ 0 on  I. 
This  condition  will  be  removed  in  Section  5. 

4»1  Representations  of  ^(a) 

3 

Recall  that  under  our  assunptims  about  a,  b,  f we  have  € C (I) 
and  hence  P - au^  < Al). 

Lenma  4,1:  Suppose  that  u^(x)  0 for  edl  x < I and  set  “ 

(x^+x._j^)/2,  ‘ idiere  |o(5.)l  = iiiax{|p(x)|  ,x<Ij}. 

Then 


(4.1) 


(1+0(K(A))),  as  K(A)  0 


» 


vA^ere  the  constant  in  the  bound  of  the  0->tem  depends  on  a,  b and  jlifiiQ  3^* 
ftoof ! Set 


(4.2)  • Pj»  X € j s l,.,.,m 

and  define  jj'Po  4 ^ E^d.)  as  the  solutions  of 
A|J  *»J  w j 

(4.3)  r ^ " if...»j'‘ 

= r,  (V).  V V « Ej<lj)  J 

respectively.  Then  we  have  ®.  = 

By  assumption  there  exists  a constant  > 0 such  that  | p(x)  I i Pq 
for  all  X < I,  and  hence  that  for  edl  sufficiently  small  fi 
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(4.4) 


!pJ 

|a.(x)|  < h.,  V X € I.  . 

J f*Q  1 3 


Mote  also  ttat  for  all  x € I. 


(4.5)  a.  s min{  a(x)  x€I-}  = as  -►  0 

Since  for  any  v i Eq(I.) 

j llvil^  -if  a(x)v*  (x)^dx](l+0(h?)),  as  h.  -►  0 , 
-j  I.  3D 


it  follows  from  (4.3)  and  (4.5)  that 


2 

•'’'^ii"E*r  I nvT" '^0 

■j  ^ Lv^Egdj)  ^j'''  E 


-2 

p. 


a. 

“3  L 


r =“p  . 

fl  J vdxl^/  /(v')^dx\li 

L ij  JJ 

[ / Cv' )^dx]  (l+0(h.)) 
-i  I.  ^ 


where  ^ solution  of 


-V”  = 1,  v(x.  ^)  = v(x.)  = 0 . 

J J 


This  implies  that 
(4.6) 


0? 


I.K.j’^E  ■ I7a7~  as  hj  - 0,  j s l,...,m  . 


■)  -.V  - — ^j.i/2 


In  order  to  estinate  cp^  j we  proceed  as  in  the  proof  of  (3.14),  The 
smallest  eigenvalue  of  the  operator  L on  is  bounded  below  by 
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a,jii^/h?.  Hence  by  (4,3)  it  follows  that 


h? 


v^iich  together  with  (4,4)  implies  that 

5? 


(4,8)  , K .!!?  s — J—  hhl+0(h.)),  as  h.  0,  j = 1,,,,^ 

Ij  2,3  E 3 3 3 


Combining  this  with  (4.6)  we  obtain— with  some  |a|  s 1- 


(4.9) 


l|(P  II 2 s l|(p  .11 2 +2a  lltp  .11  Il(|>  .1!  + II®  .11* 

I.'3E  I.'^l.j'E  I.”^1,3EI.^2,3'e  I.  2,3  E 

3 3 3 3 3 

-2 

1 3 

“ hj(1^0(hj>)),  as  h«  0,  3 “ l,,»,,in  , 

®j-i/2  3 3 3 


By  definition  (3,17)  of  0(A)  this  proves  the  lemma, 
A partition  A shall  be  a (C,m)-partition  if 


(4.10) 


^ ♦ 

5(Xj)  “ ^ , 3 " 0,1,,,,, m 


for  seme  function 


(4. 


11)  ? ( P^(I,A),  C'(x)  > 6 > 0,  V X ( Ij,  j = l,...,m,  C(0)  = 0,  C(l)  = 1 . 


Note  that  any  partition  A is  a (K,m) -partition  for  the  piecewise  linear 
function  5 ^ S(I,A)  defined  by 


(4.12) 


5(x)  “ TIP"  ^ ^ ^ ) * l,.,.,m  , 


For  a (5  ,m) -partition  we  have 
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(4.13) 


C'(t)dt  = h.^'Cx.  ,,,)(l+0(h.)),  as  h.  - 0 
3 ]-l/2  3 3 


where  the  constant  in  the  bound  of  the  0-term  depends  on  5 but  not  on  m. 

In  terms  of  (?,m) -partitions  our  Laima  4,1  can  be  rewritten  as  follows: 
Theorem  4.2:  For  the  (£,m) -partition  h we  have 


(4.14) 


(1+0(R(A))),  as  h(A)-^  0 


vdiere  the  constant  in  the  bound  of  the  0-term  depeixls  on  a,  b,  f and  c 
but  not  on  m. 

Because  p i C^(I)  and  1/5'  € T^(I,A),  the  proof  follows  directly  from 
the  fact  that  the  expression  for  Q(A)  in  Lemma  4,1  is  a Rieraann  sum  of 
(4.14). 

By  combining  Theorems  3.1,  3.2,  and  4,2  we  obtain  the  following  result: 
Theorem  4.3:  For  the  (5,m)-partition  a the  error  satisfies 


(4.15) 


(1+0(R(A))),  as  R(A)  0 


where  the  constant  in  the  bound  of  the  0-term  depends  on  a,  b,  f and  5 
but  not  on  m. 


4,2  Optimal  Partitions 

The  error  formula  of  Theorem  4,3  suggests  that  we  consider  minimizing 
the  variational  integral 


(4.16a) 
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subject  to  the  boundary  conditions 


(4.16b) 


5(0)  = 0,  5(1)  = 1 


The  Euler  equation  is  directly  solvable  and  the  functions 


(4.17) 


5(x,y)  = y / 
0 


pit) 

sSt 


2l 


|l/3 


dt,  X 5 I, 


fom  a field  of  extremals  in0<x,  y<1*  A standard  application  of  the 
E-function  test  (see,  e.g.,  Cl])  then  proves  the  foUovdng  result: 

Theorem  4.4:  For  all  functions  (4.10)  we  have 


(4.18) 


J(5)  2 J(Cg)  = 1/Yg 


v^ere  Cg  = 5(.,Yg)  with 

(4.19)  = 


0(t) 

a(t) 


2i 


1/3 


dt 


-1 


Moreover,  equality  holds  in  (4.18)  exactly  when  5 = 5g. 

Note  that  the  function  5g  belongs  to  the  class  of  functions  (4.11) i 
in  fact,  we  Vave  Cg  ^ C^(I)  and  ^gCx)  5 6 > 0 for  x € because 
lo(x)l  - Pg  in  that  interval.  For  the  partition  Ag  given  by  5g  we 
obtcun  from  Theorems  4,3  and  4.4  the  following  error  formula. 

Theorem  4.S:  The  (Cg,m) -partition  Ag  is  asymptotically  optimal  with 


(4.20) 


i”®"e 


12m‘Y 


2'  (1+0(R(Aq))),  as  R(A)  0 


where  the  constant  in  the  bound  of  the  0-term  depends  on  5g  tut  not  cxi  m. 
By  ’’asymptotic  optimality"  we  mean  here  that  for  any  other  partition 
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with  sufficiently  small  h the  error  is  larger  than  (4,20). 
For  any  (C,m)-partition  A set 


(4.21) 


1^  I.  ^ " ^»***»’" 


These  are  related  to  the  functions  of  (3.13)  ty 


(4.22)  “ ^j(^»iT')^(i‘*^(h.)),  as  hj  0,  j = l,.,,,m  . 


This  follows  directly  from  the  fact  that  the  expression  (4.9)  for  the  norms 
of  «.  is~up  to  a factor  (l+0(h.))— a Riemann  sum  of  (4.22). 

J <J 

For  the  optimal  partition  Ap  we  have 
2 1 . 


3^  ’ ■ IS?  J 


(4.23) 


■yj  / ^Q('t)dt  - g' '5'  I j ~ l»...»m  . 


ITm^Yj  I. 


12m''Y: 


In  other  words,  for  5p  all  are  exactly  equal  and— by  (4,22)— all 
j 3re  asymptotically  equal. 

ij  3 t 

Since  the  >5^  are  rot  readily  computable,  we  turn  now  to  the  quantities 


(4.24)  U.(5,m)  = _ llz.|l  , j s l,.,,,m 

j J-j  J ^ 

vMch  can  be  calculated.  For  the  optimal  partition  A^  we  obtain  from 
(3.16),  (4.13)  and  Theorem  4,5  that 
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(4.25)  j ll'V.II^  = j-g.  (l+0(h(AQ)))(0(h(Ap)),  as  R(Aq)  - o 

j ^ 12m  yq 

i^iere  the  constant  in  the  0-term  depends  r,n  a,  b,  f and  Kq  but  not  on  m. 
Thus  it  follows  from  (4,22),  (4.23)  and  (4,25)  tlat 

4.(C.,m)^  = |l<P.+t.|lp  = — (1+OC^)) 

^ ° Ij  3 3 E i2m>J 

q'/w  '5/11  y/'y  ’ If  (l+0(ll))0(h^^^) 

+ — ^ (l+0(R))0(h)  = —Vt  (1+0(R^''^))  . 

12m>J  12m'’Yj 

We  summarize  this  in  the  following  form; 

Theorem  4,6;  For  the  optimal  partition  we  have 

(4.26)  = “Vt  (l+0(R(Ag)^^^)),  as  h(Ag)  0 , 

12m  Yq 

where  the  constant  in  the  bound  of  the  0-term  depends  on  5q  but  not  on  m. 

Thus,  we  see  that  eilso  the  quantities  Uj(5Q,m)  are  a8ymptotic^dly 
equal. 

From  (4.13)  we  find  for  the  steps  hj  of  the  optimal  partition  Aq 

that 

(4.27)  h.  = (l+0(i)),  as  m - 


and  hence  that 
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hCAg) 


> \q  > 


0 


with  a constant  Xq  tlmt  does  not  depend  on  n but  only  on  the  problem 
and  Pq.  In  other  words,  the  optijnal  partition  is  (Xq,1) -regular. 

Conversely,  it  turns  out  that  asymptotically  the  optiiral  partition 
Aq  is  characterized  by  the  asymptotic  equality  of  the  |j,j.  This  is  the 
content  of  the  following  theorem. 

Theorem  4.7:  For  the  partition  A suppose  that 

(4.28)  ||z.(A)||^  = n(l+0(h^^^)),  as  h(A)  ^0,  j = l,...,m, 

■^3  ^ 

where  p does  not  depend  on  j.  Then 


(4.29)  j||e(A)||j,  = j||e(AQ)iyi+0(h-''^)),  as  h(A)  --  0,  j = l,...,m, 
and 

(4.30)  |Xj-Xj^|  = 0(h^^^),  as  h(A)  -+  0,  j = l,...,m. 

Proof:  We  show  first  that  A is  (X,l)-regular.  For  ease  of  notation  let 


- [i  "3  1 


1/2 


, j - 1,... ,m. 

For  any  j = l,...,m,  we  have  by  (4.9),  with  some  |a|  s 1, 


(4.31) 


Tl?h?(l*0(h))  + 2ar,.h?'2  jJI+jllj.(l+0(R))  t J. 
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Let  now  h.  = h.  Then  (3.16)  and  (3.6)  show  that 
^0 

I.  If ll"llE==h2 


30  ^0 


'0 


or 


T llV-illp 


- >;3/2  r . 

- h 0(h  ),  as  h -*•  0 . 


(4.32) 

Hence  we  obtain  from  (4.31)  and  (4.28)  that 


or 


(4.33) 


4^(l+0(h^''^))=  T)?  h^(l+0(h))  + h^  0(h^^^)  + h^  0(h) 

3o 


= Ti?  R^(l+0(h^^^)),  as  h ->•  0 . 
^0 


Now  let  h.  = h.  Then  we  have  instead  of  (4.32) 


3 


j ||^|/  11^  = h h^/^0(h^^^)  , 

3l  1 


and  hence  by  (4.31),  (4.28),  and  (4.33) 
2 r3/, ,„,rl/2^^  _ 2 ^3, 


5/2  rl/2Q^p^l/2j  ^ j^2pj 


(4.34)  nf  r(l+0(r'^))  = Tif  h^(l+0(h))  + 2crn.  h h' 

3o  h h 


In  other  words,  z = h/h  satisfies  the  polynomial  equation 
(4.35)  q?  (l+0(h^^^))z^  + 0(R)z^  + 0(R^''^)z  - q?  (l+0(h))  = 0 . 

3o  h 

By  comparing  (4.35)  with  the  equation 


(4.36) 


q?  (l+0(h^^^))z®  - q?  (l+0(h))  = 0 , 

3o 


Rouc’ie's  theorem  shows  that  for  sufficiently  snail  h we  have 
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(4.37) 


This  shows  that  A is  indeed  (X,i)-r^lilar. 

Our  assunqptions  about  a and  p ensure  tiiat  for  any  fixed 


3, 


We  represent 
Theorems  3.1, 

(4.38) 


and  hence,  by 
(4.39) 


A by  the  piecewise  linear  function  ? 
3.2,  and  Lemma  4.1  it  follows  that 

jlie|||  = Ti^h|j  (l+0(h))  ^ O^h^  , 

(3.16)  that 

j < Qu  j j|el||  < Oi^hjh^  5 (h^hjh 

^ D 


of  (4.12). 


» 


Tiien  fixan 


where  in  the  last  inequality  the  regularity  of  A was  used.  Therefore,  we 
have 


jUJL  = Ti.hy^0(h^^^) 

j J J J 

whence  by  (4.31) 

I llzJlp  = Ti?h?[(l+0(h))+0(h^^^)+0(h)]  = Ti^h?(l+0(h^^^))  , 

j J ^ J J J J 

that  is,  by  (4.28) 

= qjh?(l+0(h^^^))  . 


(4.40) 
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Since 


2 1/3 

) dt|  (1-KKh)),  as  h -»  G , 


(4.40)  ia;>lies  that 


(4.41) 

which  by  (4.19)  shous  Hat 

(4.42) 


dt  = (12ii*)^^(l40(h^^)) 


ro  0 L^J 


* 2 -1/2 
dt  = m(12n'^)  (l+0(h^'^)) 


By  (4.38)  and  (4.uo)  we  obtain  now 

jlleli^  = n^i^(l+0(h^^^))  , 

vduch  by  (4.42)  and  Theoreni  4.5  gives  (4.29).  Finally  ftora  (4.41)  it  follows 
sunnation  over  the  first  j intervals 


rrl/2. 


/ * I I <it  = zrr-  (l+0(h"' "))  . 


l_atTj 


ntTr 


On  the  other  hand,  we  obtain  from  (4.23) 

^0 


Thus  we  have 


/ 

x40 
: 

which  duplies  (4.30), 


[^r 


mrr 


= 0(h^^^)  , 
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Sujjpose  that  the  optimal  partition  function  is  changed  to  some 
partiticai  function  5 + e satisfying  (4.11).  By  (4.13a)  the  derivatives 

control  the  stepsizes,  and  hence  we  assume  that  is  small.  For  any 

given  m,  let  jlkllg  ^nd  errors  associated  with  the  (?,m)- 

partiticai  and  (5Q,m)-partition,  respectively.  Then  by  Theorem  4.3  and 
(4.16a)  we  have 


|j!|e|ij,  - I ~ » as  m . 


Since  the  variational  int^ral  J is  stationary  at  we  have  = 0 

and  hence  it  follows  from  the  meaiv-value  theorem  that 

|j(?)-J(4o)|  = 0(j|leM!V 

Therefore,  because  (4.15)  inplies  that 

lllell^  + jllegllj,  = as  m - - , 

we  obtain 

Ijllellj;  - iII®oIIeI  ■ I**^'*^c  ^ * m 

In  other  words,  a change  of  Cg  ly  some  small  leads  to  a change  of 

2 

the  error  proportional  to  ^||e'||^.  This  shows  that  the  value  of  the  optimal 
error  is  rather  stable  ui.der  perturbations  of  the  optimal  partition.  On  the 
other  hand,  the  optinal  partition  itself  is  not  too  stable  and  hence  needs 

to  be  computed  only  with  relatively  low  accuracy.  ; 

\ 

By  Theorem  4.7  the  optinal  partition  is  clarHcterized  by  the  asymptotic 
equality  of  the  quantities  = |i.j(5g,m)  of  (4.24).  Let  r^  agair.  denote 
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the  residual  r = lXu.)-f  (m  the  subinterval  I.  airKi  set 

o 3 

C‘»-®»3)  w|  = / r-(x)^dx,  j = 1,...^  . 

Ve  call  the  quantities 
2 1 

rt.W  5 = 1.2,...^ 

the  error  indicatxrs  fo*.  the  intervals  L , . . . ,I  , and  set 

1 IS 

c 

(^.45)  £(A)  = { £ e . 

3=1  ^ 

All  these  quantities  are  directly  ooeputahle  once  the  finite  eleaesit  solutioR 
is  know.  Ihe  next  theorem  shows  that  e(A)  is  asyEiptotically  equal  to  the 

IheoreB  4.8:  a)  For  any  iJartition  A ve  have 
(4.46a)  jie|£  = e(A)^(l*0(h)),  as  h(A)  - 0 . 

b)  If  A is  (X,»c)-regular,  1 ^ k < 2,  then 
(4.46b)  V?  = p|h-(l+0(h'^)),  as  h(A)  -►  0,  3 = , 

tAere  6 = 1-  ic/2. 

Proof:  a)  The  definitim  of  p-  given  in  hsaaa  4.1  provides  that 

jp.j  = rax  }o-(x)|  > p.  ^ 0 , 

^ x€l.  ^ 

3 

and  by  (4.4),  (3.15),  and  (3.6)  we  have 
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Hence  by  (3.11)  and  (4.2)  we  obtain  with  some  ja^j  5 1,  i = 1,2,3,  that 
V?  = / (p.+a.(x)+T-(x))^dx 

n * *1  T 


(4.47) 


3 x'  3 
"3 


= . i 5|hjo(h ) . j iit  nj 

^ Pn  3-' 


0 

* ^®1  ^ ^ 2a2Pjhy^  ^“3  ^ 


With 


-2.  3 

A 4 m * o o 

S(4)2  = A I -J-L  . T(4)2  = Z IIt-II^ 

12  j=i  a._^2  j=l  Ij  3 u 


it  then  follows  tlat 


e(A)^  = S(A)2  ♦ ^ S(A)2o(h^)  + T(A)^0(h^) 


+ 2 ^ S(A)^0(h)  ♦ 2(i2S(A)  T(A)  0(h)  ♦ 2a3S(A)T(A)0(F’) 
This  proves  (4.46a)  since  by  lama  4.1  and  Iheorens  3.1  and  3.2 


S(A)^  = TieSr(l+0(h)),  T(A)^  ^ c I T = C -rleip  . 


b)  For  (X,K)-regular  A with  1 5 x < 2 we  have  by  (3.6) 


I J'fjSo  - ^ I J^E  - ^ ^ b^ 


Hence  (4.47)  igplies  tJat 


i = pji.[l+  -4  O(n^)  ♦ \ * ~ 0(h)  ♦ — OCh*')  + -4  0(h^^)  j 


'j  - 7 

^0 


2 

^0 


which  is  (4.46b)- 
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It  nay  be  noted  that  in  [3]  we  proved  an  iq>per  bound  for  jllellg  of 

2 

the  form  (4.46b)  with  1/12  replaced  by  the  larger  factor  l/n  . 

Theoran  4.8  states  that 

lllellr 

(4.48)  0(A)  = = 1 + 0(h(A)>,  as  h(A)  - 0 . 

In  other  words,  the  effectivity  quotient  0 tends  to  one  with  R ->  0. 

2 2 

In  contrast,  the  corresponding  estimates  in  [3]  only  provided  for  0 < 12/n  . 

We  e3q)ect  the  error  indicators  to  be  asyirptotically  eqtal  to  the 
quantities  of  (4.24).  Theoran  4.9  below  shows  that  this  is  indeed 
correct  for  regular  partitions.  Hence  our  aim  is  to  construct  such 
partitions  for  which  all  6j  are  asymptotically  equal.  It  turns  out  that 
—as  before— these  partitions  are  automatically  regular. 

Theorem  4.9;  (a)  Let  A be  a (X,ic)-regular  partition  with  1 s k < 2. 

Then 

(4.49)  T ||z.||p  = 6?(l+0(h^)),  as  h(A)  0 , 

■•■j  j ^ ] 

with  e = 1 - k/2. 

(b)  All  partitions  A for  which 

(4.50)  6j  = M-Cl+od))  as  h(A)  0,  j = l,...,m  , 
with  n independent  of  j,  are  (X,l) -regular. 

Proof;  As  in  the  case  of  (4.39) it  follows  from  (3.16)  and  the  assumed 
regularity  of  A that 

II.  ii2  ^ « 2,2r2  , „ 2,3r25 

I JI+jllE  = C r,jh.h  < C n^h-h  , 


and  thus 
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T ll’l'-;llp  = -n-h^^^OCh^),  as  h(A)  0 . 

-■■j  J ^ 3 ] 

Now  (4.31)  shows  that 

\\z.\\l  = Ti^h^d+OCh))  + 2oniMo(h^)(l+0(h))  + Ti?h!o(h^^) 

33  3 3 

= Ti?h?(l+0(h^))  . 

Since  by  Theorem  4.8(b) 

-Z3 

(4.51)  Ti?h!  = A p ^ = ej(:+0(h^))  , 

this  proves  (4.49). 

(b)  Because,  generally. 


it  follows  from  (4.47)  that 


'3  ^ * Oih^^h 


whence 


a a*  0(1))  = Cj  = + h?IiO(h^^)  . 

Now  si5)pose  that  (4.50)  holds  for  t and  tlat  h.  = R;  then 

3o 

H^(l+o(l))  = h^(Ti?  +0(h^^))  . 

^0 

Similarly,  for  h.  = h we  obtain 
d " 

ii^a*  0(1))  = h?  h^  + h^(h^^^) 

3-!  ~ 


- 28  - 


and  hence  z = R/h  satisfies 

[t)?  +0(R^'^^)3z^  = ti?  + 0(h^^^)z  . 

By  ccjtiparing  this  with  the  polynomial  for  z in  which  the  last  term  on  the 
right  has  been  dropped,  it  follows  by  Rouche's  theorem  that 


for  sufficiently  snail  h and  therefore  that  A is  (X,l) -regular. 

Theorems  4.9  and  4.7  confirm  that,  as  expected,  our  aim  should  be  to 
construct  partitions  for  which  all  Cj  are  asyirptotically  equal.  Then  the 
error  of  the  partition  will  be  close  to  the  asyirptotically  optimal  error 
(4.20).  A natural  approach  for  this  construction  is  the  use  of  an  adaptive 
mesh  refinement  algorithn  of  the  form  discussed,  for  instance,  in  [4].  We 
shall  not  repeat  the  details. 
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5.  Optimal  Partitions  - Part  II 

In  the  previous  section  we  assumed  that  Ug(x)  i 0 for  x € I.  dearly, 
this  represents  a severe  restriction.  Actually,  the  results  are  largely 
valid  also  vdien  uJJ  has  zercfjs  in  I,  but  the  proofs  becone  more  delicate. 

We  illustrate  the  approach  for  the  frequent  case  vdien  uJJ  € C^(I)  has  finitely 
nany  simple  roots  in  I,  say, 

(5.1)  ujcy  = 0,  u'”(y  0,  k = l,...,q,  0<  ^<  ^<  ...  < ^<1. 

Lemma  5.1;  (Aider  the  stated  cissvnrptions  we  have  for  any  (X,K)-regular  parti- 
tion A with  1 < K < 2, 

(5.2)  Q(A)^  = A I h?  (1+0(R(A)®)),  as  ii(A)  0 

^j-l/2  ^ J • 

vdiere  e = 1 - ic/2  and  the  constant  in  the  bound  of  the  O-terni  depends  on 
a,b  and  f . 

Proof:  Because  of  (5.1),  we  may  choose  C2>>  > 0,  6q  > 0,  such  that 

(5.3)  C2|x-Cj^|  > |p(x)|  > V |x-^|  < 6q,  k = l,...,q  . 

For  any  6 > 0 we  introduce  the  sets 

Ig  = {xa  |x-^|  < 6 for  some  » 1°  = IS  Ig  , 

Jg  = {j=l I*  n ^ m}  n Jg  . 


(5.4) 
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-1 


We  assume  that  6q  2 (8q)  and  hence  that 


hj  < 2(6Q+h)q  < 46Qq  < y , for  h < 6^  . 


Since 


min  { I p(x)  1 ,x6lc  } = > 0 , 

0 

the  estimte  (4.9)  of  the  proof  of  Lemma  4.1  holds  for  the  subintervals 

I.  with  j € J?  ; tMt  is 
^ 0 

^ llcp.d  = Ti-h^(l+0(h))  as  h(A)  0,  j € . 

J-j  3 4-  ] 3 Oq 


Hence  for  h S 6q  we  have 


(5.5)  Q(A)^  ^ t I ^ Ti2h^](l+0(h))  > Cp^xV^  I n 

3-1  u . ,c  3 


2^2r2< 


■rr2K 


jO 


j€J 


ric/2 


Now  consider  the  sets  (5.4)  with  6 = h 5 6^.  Then  (5.3)  implies  that 


|S.|  > J JC  _ 


and  hence  (4.8)  modifies  to 


(5.6)  T ||(p„  .\\l  < C iTiM(l+0(h))  < Cn?h?h^^  as  h(A)  - 0,  j € . 

^ J J J J 0 


On  the  other  hand  we  have 


|a-(x)|  5 2 nex  |p(x)|  5 Ch'^^^,  j € Jc 

^ x€I.  ° 

3 


whence  by  (4.7)  and  (5.5) 
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(5.7)  T ||(p,  aI  S Ch?h'^(l+0(h))  < CQ(A)¥%,  as  h(A)  0,  j € J-  . 
J-j  J J ® 

For  ease  of  notation  wet 


r\  ni  O 0 0 ^ O 

S(4)2  = .1  l^h?.  Mi)  = Z jjk  .||2  . 

j=l  J •>  T=1  j 


3=1  " 3=1  3 

Then  it  follows  from  (5.6)  and  (5.7)  that 


(5.8) 


R(A)^  < C[  I + 0Q(A)^h^^  I h. 

3€J,  ^ 


< C[S(A)^+Q(A)^]h^® 


In  the  case  of  Lenma  4.1  the  assumption  Uq  = 0 does  not  enter  into 
the  proof  of  (4.6),  and  thus  we  have  also  in  the  present  case 


(5.9)  J ||(P^  -lip  = Ti?h?(l+0(h)),  as  h(A)  -►  0,  j = l,...,m 
Ij  1>J  E 3 3 

and  therefore 

lU  ^ 

(5.10)  y T K = S(A)''(l40(h)),  as  h - 0 . 

j=l  E 

Altogether,  with  some  suitable  constants  o,B  € [-3,1]  it  follows  fixsn 
(5.8)  and  (5.10)  that 


o ni  - / m ,'.2/m 

= S(A)^(1+0(S))  ♦ aCS(A)(l+0(H))(S(A)^+Q(A)^)^^' 


.1/2 


♦ PC(S(A)^+Q(A)^)h^®  . 


R(A)^ 


After  separating  the  middle  term  or.  the  right  and  squaring,  we  obtain  the 
equation 
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CS(A)^(S(A)^+Q(/^)W® 

= Q(A)'*(l+0(h^®))-20(A)^S(A)^(l+0(H^®))  + S(A)^(1+0(R^®))  , 
which  has  the  solutions 

Q(A)^  = S(A)^(1+0(H^®))  ± [S(A)VR^®)]^^^,  as  h(A)  - 0 , 


and  hence  proves  (5.2). 

Now  the  theory  of  Section  4 can  be  carried  over  to  the 
present  case.  As  before,  we  consider  partition  functions  C,  but  here  we 
need  to  weaken  (4.11)  by  requiring  instead  of  C*(x)  2:  6 > 0 on  each  I. 
t)«t(p(x)/C'(x))  isRiemann  integrable  on  each  subinterval.  Moreover,  we 
assume  that  for  given  K and  n • the  resulting  (C  ,m)-partitions 

are  (X,K)-regular  with  1 s ic  < 2. 

Then  as  in  the  case  of  Theorem  4.2  it  follows  t^at 


(5.11) 


/ p(x)  Y 1 

I a?xT 


(1+0(H(A)®)),  as  h(A)  - 0 , 


vihere  the  constant  in  the  bound  of  the  0-tem  depeivds  on  a,  b,  f and  5 but 
not  on  m.  This  suggests  again  oonsideration  of  the  variational  problems 
(4.16a/b)  and  hence  of  the  optinal  partition  function  of  Theorem  4.4. 
□.early 


g^'xT  ~ ^ (p(x)a(x))^'^ 

^0'  ^ “0 


o 

o 


is  Riemann->integrable  on  f. 
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We  show  first  that  the  (5Q,m)-partitions  are  (X, 5/3) -regular.  Since 

(5.12)  lp(x)|  > C min  {l,ndn|x-^,  | },  V x € I , 

k ^ 

with  some  C > 0,  it  follows  that 


Cihj  > / 5j(x)c3x  = ^ ^ C2  / p(x)^^^dx  > Cghj^^ 

I • I ♦ 


- -3/5 

where  again  > 0.  This  implies  that  h 2 Qm  as  well  as  1/m  5 ()h,  and 

hence  the  partition  is  indeed  (X,5/3)-regular.  Now  Theorem  4.5  is  easily 

— — 1/G 

shown  to  hold  with  0(h)  replaced  by  0(h  ). 


r5/6 


For  any  which  intersects  Ig  with  6 = h we  can  use  (5.6) 
in  the  estinates  leading  to  Theorem  to  obtain 


(5.13)  as  h ^ 0,  j e 6 = h^^®  . 

^ 12m  Yg 

In  other  words,  for  the  optimal  partition  the  are  asymptotically  equal 
for  all  intervals  which  are  not  too  close  to  a root  of  Ug.  As  the  numerical 
examples  of  Section  6 show,  the  ^ij(^Q,m)  for  the  intervals  close  to  roots 
are  generally  larger  ajid  the  ratio  of  the  largest  to  the  smallest  of  these 
values  does  not  tend  to  one  for  m 

The  analog  of  Th  .7  is  somewhat  more  complicated.  We  formulate 

it  as  the  following  theorem. 

Theorem  5.2;  For  the  partition  A suppose  that 

(5.14)  j |{z.(a)||j,  = |j.(l+0(h^)),  as  h(A)  0,  j = l,...,m  , 

j ^ 

where  (i  does  not  depend  on  j and  e = 1/12.  Then  A is  (X, 5/ 3) -regular  and 

(5.15)  j||e(A)||£  i jlie(Ag)il£  (l+0(h^)),  as  h(A)  -*•  0 . 
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Proof:  1)  We  show  first  that 


(5.16) 


j — Ch.  , 

j J J 


with  sane  C > 0 which  depends  only  on  p.  From  (3.13)  it  follows  that 

iJkjllj;  = sup/— I / p(x)v(x)dx|  V V € Eq(I.)|.. 
j Llj  E I.  J 

Clearly,  for  given  y,6  the  function 

Vj(x)  = Yhj  + ^ 6h?)(x-Xj_j^)  - Y Y(x-Xj_^)^  - i 6(x-Xj_^)^,  V x € 

belongs  to  EQ(Ij)  and  a short  calculation  shows  that 

C,a?h?  < / v'(x)^dx  < C.a?h?,  / v(x)^dx  S C.a?h! 

1 3 3 T 2 3 3’“’  333 


where 


a.  = (r^+(6h.)^)^^^ 


and  all  constants  are  positive.  Hence  also 


Now  with 


Y = p(Xj_^),  6 = Pj(Xj_j^),  p(x)  = -vl!(x)  + Pj(x),  V X € 


|p.(x)|  = o(h.),  as  h.  -»■  0 , 
J J D 


we  obtain 
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I / p(x)v(x)dx|  > / v'(x)^dx  - ( / p-(x)^dx)^^^( 

I.  I.  I-  ^ 

D 3 3 

v^ence — for  sufficiently  srifill  h — 

3 3 

(5.17) 


By  assuirption  we  have 

(5.18)  u5(x)^+  uj''(x)^>  C > 0,  V X € I 
and  hence  also  with  some  C > 0 

“ h?(p(Xj_^)^+p’(Xj_^)^)  ^ C3i?  > 

which  together  with  (5.17)  proves  (5.16). 

2)  Next  we  show  that  for  sufficiently  small  hj 

(5.19)  jllell^  > j llellj,  > C h^^^  C>0. 

Obviously,  we  have 

T ||e||p  > inf  / a(x)[ui(x)-v'(x)]^dx 
j V I.  ^ 

where  the  infiraum  is  taken  over  all  linear  functions  on  I. 


/ v(x)^dx)^^ 


Thus  v'(x) 


is  constant  and  it  follows  that 
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, mI  ^ CC  / ul(x)^dx  - ^ / u;(x)dx)^] 


‘j  ^ 


hj  “0- 


Since 


uj(x)  = uJ(xj.^)  ♦ u5(x..;i)(x-X5.i)  ♦ |uJ”(Xj.j)(x-x.)' 


+ o(hf),  as  h . -►  0 , 


a single  ccLLculatLon  shows  that 


j ||e|||  > C[u5(x..^)^h?  ♦ uj"{x..^)\5]  - o(l)h|  . 

Using  (5.18)  we  then  obtain  (5.19)  for  sufficiently  small  h^. 

3)  Now  we  can  show  that  A is  a (\,5/3)~rsgular  partition.  By  TheoresD  3.1 
and  (5.14)  it  follows  that 


(5.20) 


jllell^  = n*i^(l+(h^'^^))  , 


and  thus  b'/  (5.19)  that 


(5.21) 


r ^ „ 1/5  2/5 
h < C m p. 


On  the  other  hand,  we  have  by  (5.9),  (4.7),  (4.4),  and  (3.16) 


p(l+0(h^^^^))  = 


(5.22) 


5 CC-n-h^'^^+h^'^^+h.  j ||e|L]  < C[h?^^+m^'^%.] 


Let  h.  = h and  note  that  m i 1/h,  then 

3o  - 
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and  hence  for  sufficiently  small  h 


(5.23) 


h > C 4 


2/3 


On  tile  other  hand,  (5.16)  and  (3.16)  give 


H(l+0(h^^^^))  > ^ i.li'PjllE  - 

> C(hf  2-h.  . !le|L) 

3 3 -j  ^ 

and  tiierefore,  by  (5.21)  and  (5.23),  for  h.  = h 

3o 


(5.24) 


H > > cCh^' 


If  the  term  on  the  right  is  positive  then,  with  some  positive  constant 

—1/2 

(independent  of  h),  z = h must  satisfy 


5 2/3  2 _ 

z - n z - Cji  . 


^ Cauchy’s  rule  this  iii?>lies  that 

since  by  (5.23)  we  have  n ^ 1 for  sufficiently  small  h.  On  the  ocher  hand, 

if  the  tern  on  the  right  of  (5.24)  is  negative,  then  we  have  inmediately 
4/9  2/5 

h ^ (i  < [i  . Thus  with  (5.23)  this  gives  indeed 

h^^^  < C h . 


•*)  For  the  proof  of  (5.15)  let  now  6 = h . By  (5.12)  we  have  > C6  for 
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j € Jg  anu  hence 


JA  ■ '(?)’■  > 


Thus  from  Lenna  5.1  and  Theorems  3.1,3. 2 it  follows  that 


(5.25) 

and  therefore  (3.16)  gives 

(5.26) 


j||e||2  i C i n-h^  » 3 « . 


iPik  ^ c h.K  < c ^-1  f 3 


respectively.  Note  that  in  the  second  inequality  the  (\, 5/ 3) -regularity 
of  A vas  used.  Together  with  (5.6)  and  (5.9),  (5,26)  leads  to 

and  therefore  by  (5.14)  to 


(5.27) 


= Ti?h?(l+0(h®)),  j ^ Jg  . 


Hence  analogous  to  (4.41)  we  have 


(5.28)  / = (12H^)^^^(l+0(h®)),  j € Jg  . 

Let  and  m2  denote  the  cardinalities  of  Jg  and  Jg,  respec- 
tively, that  is,  m = + m2.  We  want  to  show  that 

"^2  -e 

— = 0(6;  = O(h^)  . 

mi 


(5.29) 
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For  this  note  first  that  because  of  < C6  for  j € Jg  we  have  by  (5.22) 


(5.30) 


. 5/2^,  1/2  . . , , 

ji<C(6h.  +h.  +h.m  n),  j € Jc 
3 3 3 ’ 6 


and  thus 

(5.31) 

Suppose  first  that 


h.  > C min 
3 


j f J,. 


2/3 


then 


(5.32) 


.1 

nu  < 6 (min  h.)  < C ^557*-  , 3 € J-. 

^ j€J  3 ^2/3  6 


Because  of  > (36  for  j € Jg,  it  follows  from  (5.27)  that 


/ \2/3 

V<=(t)  • 


3 ^ h 


and  hence 


-1  / fi  \ 2/3 

nu  > (l-6)(maxh.)  > c(-) 


Ttogether  with  (5.32)  this  proves  (5.29)  in  this  case. 

2/5 

Now  suppose  that  in  (5.31)  we  have  hj  > (3^,  . Then 


m2  ^ C!6|i 


-2/5 


while  (5.24)  ijiplies  that 
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m ^ ”2/5 

> Cji 


Together  these  estinates  show  once  more  -Oiat  (5.29)  holds. 


Finally  consider  the  case 


h.  > j € Jg. 


1/9  1/9 

in2  < C6m-^  = CSCm^-hn^) 


n^-  C%^m2  S . 


By  applying  Rouche's  theorem  to  the  pair  of  polynomials 


2 2 2 2 
z - az  - = 0,  z - am^  = 0,  a = C 6 

it  follows  readily  that 


1/2 

5 C6n^'  2 C6nij^ 


Thus  (5.29)  is  valid  and  we  have 


(5.33) 


= m(l+0(h^''^^)) 


By  definition  of  it  follows  now  from  (5.28)  and  (5.33)  that 

0 


, 1 r ,^x2t1/3  r 1/3 

^ ^ i [^  ] ^ ik  ij  N ] '' ' h i 


r 7 "I 

r P(t)^ 

^ iTtT 

• L.  • 


2T1/3 


> ni^(12u^)^^^(l+0(h^)  = m(12n^)^''^(l+0(h®)) 


i 
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or 


(i+o(h^))  . 

Ygin 

From  (5.20)  and  Theorem  4.5  (modified  to  the  present  case)  this  implies  (5.15). 
It  nay  be  noted  that  the  analogous  relation  to  (4.30),  namely, 

[Xj-Xj  I = 0(h  ) 

is  easily  proved  when  there  is  only  one  root  of  Uq  in  1.  In  general,  the 
situation  appears  to  be  more  conplicated. 

Now  we  turn  to  the  analog  of  Theorem  4.8. 

Theorem  5.  3;  For  any  (X,k) -regular  partition  with  1 < tc  < 2 

,2.2  ^ 

(l+0(h®)) 


(5.34) 


I«"»E 


, r m vfh^ 

= -L  T _iJL 

[j=l  ^j-1/2  _ 


where  is  given  by  (4.44)  and  6=1-  k/2. 
Proof:  Recall  that 


2 2 

o T m v.h. 

e(A)2  = J^  r 


j=l  ^j-1/2 


Using 


= / (p-+a.(t)+T.(t))^dt 
J I^  J J j 

(5.35)  3 


|a.(x)|  < Ch.,  j_||t.||q  < C j jjellj. 


we  obtain,  with  certain  |C^|  5 C,  i = 1,...,5, 
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12e(A)^ 


-'^k3 

f ^ 

j=l  ^3-1/2 


h 


2 

j 


+ C 


3 


II  II  J/2 

j • 

3 


By  Lenna  5.1  and  Theorems  3.1  and  3.2  the  first  term  is  asymptotically  equal 

2 “K 

to  12  jllelij,  and  fron  (5.5)  it  follows  that  jl|e||g  2 Ch  . The  other  terms  are 
then  easily  estijiated  to  give — with  different  constants  — 


e(A)^  = j||e||^(l+0(h®))  + j||e||^  h*^^  + j||e||^  h^ 


iNIe  * C„,M^ 


,1/2 


= j||e|||(l+0(h^))  , 
which  proves  (5.34). 

Finally,  we  show  that  also  Theorem  4.9  carries  over  to  this  case. 
Theorem  5.4;  For  the  partition  A suppose  that 

(5.36)  6j  = n(l+0(h^))  as  h(A)  -*-0,  j = l,...,m  , 

where  (j,  does  not  depend  on  3 and  e = 1/12.  Then  A is  (X, 5/ 3) -regular 
and 


(5.37)  j||e  (A)llj,  5 il|e(AQ)||j,  (l+0(h®)),  as  m ->  « . 


Proof:  We  show  first  that  -A  is  (X, 5/ 3) -regular.  For  this  note  that  generally 
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j llZjllg  5 e?(l+0(hj)),  as  hj  0,  j = . 

This  follows  ftom  (3.8)  in  the  same  manner  as  (3.14)  follows  from  (3.13b). 
Thus  ty  Theorem  3.1  and  (5.36)  we  have 


jllell^  < C , 


whence  by  (5.19) 


c ^ p 1/5  2/5 
n < C m (i 


Now  (5.35)  gives 


(5.38) 


Vj  = i.lIPj+t^j+T^jllo  - C(p.hy^+h?^^+  j llellj.) 


and  thus— as  in  the  case  of  (5.23) — for  sufficiently  small  h 


h > . 


By  (5.12)  we  have 


/ p(x)^dx  > C h? 
r j 


and  tiius 


leads  to  (5.24).  The  renainder  of  the  proof  of  the  regularity  of  A now 

proceeds  exactly  as  part  3)  of  the  proof  of  Theorem  5.3. 

Similarly  the  proof  of  (5.37)  follows  that  of  (5.15).  In  fact,  for 
• c 

j € Jg  we  have 


v^iich  as  in  (4.47)  leads  to 


V?  = p?h|(l+0(R®)),  j € , 

that  is, 

= ri?h?(l+0(h®)),  j 6 . 

J J 0 

Therefore  (5.28)  holds  again,  fforeover,  because  of  < Cn^  we  obtain 
from  (5.38)  the  estimates  (5.30)  and  (5.31)  which  in  turn  inply  (5.29). 
Now  the  renaining  conclusions  of  the  proof  of  Theorem  5.4  apply  verbatim. 


- 45  - 


6.  Numerical  Exanples 

We  illustrate  the  theorerical  r-esults  with  some  cxmputatioral  results 
for  the  following  two  sample  probleas,.: 

Sample  Problem  A: 

(6.1a)  - ^ (x+a)P  ^ + (x+a)%  = f,  0 < x < 1,  a > 0, 

(6.1b)  u(0)  = u(l)  = 0 

v^ere  f is  chosen  such  that  the  solution  of  (6.1a/b)  is 

(6.1c)  " (x+a)’^  - [a^(l-x)+(l+c)^x],  x € I . 

Here  the  cx)efficient  functions  and  f ai-e  analytic  ir  I , and  we  have 
u|J(x)  i 0 for  X 6 I.  Hence  the  theory  of  Section  4 applies.  Note  that  for 
snail  a and  negative  r we  <an  create  severe  near-singularities. 

Sample  Problem  B; 

(6.2a)  -u"  + u=f,  0<x<l, 

(6.2b)  u(0)  = u(l)  = 0 

where  f is  chosen  such  that  the  true  solution  is 

(6.2c)  " e°^(x-P)  + [p(l-x)-e°(l-p)x],  a ^ 0,  p = ^ ~ . 

Here  f is  analytic  on  I and  uJJ(x)  tes  a simple  root  at  x = 1/2,  and 
hence  we  can  apply  the  theory  of  Section  5. 

The  tables  of  computational  results  given  below  include  the  following 


data: 
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m number  of  intervals  used  in  the  partition 

E = 100.|!e|L/jl|uQ(L  relative  error  in  the  energy  norm  expressed 
^ ^ in  percent 


E 


0 


1 

12m^Yo 


asymptotically  smallest  error  achieveable 
with  meshes  of  m intervals 


0 = j||e||j./6(A) 


eff  activity  quotient  (4.39) 


(tt 


( max  $.)/(  min  e.) 
j*l).«.)m  ]“l)...jm 


ratio  between  the  largest  and 
smallest  value  of  the  error  indi- 
cators 6^(A)  of  (4.38a) 


partition  points  of  the  particular  mesh 


Tables  1 through  6 concern  two  cases  of  san5)le  problem  A.  Each  time  the 

left  endpoint  x = 0 of  the  interval  is  near  a singularity  of  Uq  and  this 

is  reflected  in  the  fact  that  the  largest  and  smallest  error  indicators  6j(A) 

always  occur  on  the  first  and  last  subinterval  and  respectively. 

2 

But  in  the  second  case  the  energy  expression  includes  a weight  (x+a)  which 
goes  strongly  down  near  x = 0.  Thus  in  this  case  the  near  singularity  of 
the  solution  shows  up  more  weaJcLy  under  the  energy  norm. 

In  all  cases  the  effectivity  quotient  is  less  than  one  and  hence  the 
estimate  e(A)  turns  out  to  be  an  upper  bound  of  jllellj.*  Of  course,  the 
theory  is  only  asymptotic  in  nature  and  thus  6(A)  could  be  smaller  than 
jllellj,.  Note  that  for  relative  accuracies  better  than  10%  the  estimate  never 
overshoots  the  true  error  by  more  than  10%.  In  fact,  for  hi^er  accuracies 
0 equals  one  for  all  practical  purposes.  This  is  in  complete  agreement  with 
the  theory  and  shows  that  the  a-posteriori  error  estimate  is  very  reliable 
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and  not  at  all  pessimistic. 

In  the  presence  of  the  near  singularity,  the  use  of  nonuniform 
meshes  is  very  advantageous  (see,  e.g. , Table  6),  and  the  approximately 
optiial  meshes  produce  errors  close  to  the  optimal  values  which  by 
Theorem  4.5  decrease  with  1/m.  Here  the  "weaker”  singularity  of  the 
second  case  is  rather  noticeable.  The  nonuniform  meshes  are  only  approxi- 
mately optinal  as  the  ratio  co  shows  which  in  each  case  is  reasonably 
close  to  one  but  certainly  not  equal  to  it.  Nevertheless,  as  expected, 
the  corresponding  errors  are  clearly  not  very  sensitive  to  such  changes 
of  ‘the  mesh  except  for  low  accuracies. 

Tables  7 through  11  contain  results  for  two  cases  of  sample  problem  B. 
Essentially  all  aspects  are  the  same  as  for  problem  A.  However,  in  all 
cases  the  naximel  occurs  in  the  neighborhood  of  the  root  Xq  = 1/2 
of  u'^(x)-,  and,  as  expected,  the  ratio  co  does  not  converge  to  one.  How- 
ever, if  CO  is  conputed  only  for  all  intervals  outside  a snail  neighborhood 


of  Xg,  then  we  have  again  the  desired  convergence  of  co  to  one. 
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Table  1 


Problem  A with  p = 0,  q = 1,  r = -1/4,  a = 1/100 
Uniform  mesh,  jIIuqII^.  = 6.09811 


m 

E 

6 

CO 

5 

85.301 

22.613 

.1706 

8.84(+6) 

10 

73.768 

11.306 

.2950 

2.34(+7) 

20 

58.784 

5.653 

.4702 

5.31(+7) 

40 

41.933 

2.827 

.6708. 

l.ll(+8) 

80 

26. 3xb 

1.413 

.8419 

2.19(+8) 

Table  2 

Problem  A with  p = 0,  q = 1,  r = -1/4,  a = 1/100 
/^pioxiinately  optimal  mesh,  jliugllj.  “ 6.09811 


m 

E 

^0 

e 

CO 

5 

22.243 

22.613* 

.6524 

5.854 

10 

11.289 

11.306 

.9025 

2.274 

20 

5.652 

5.653 

.9757 

1.372 

40 

2.826 

2.827 

.9940 

1.111 

80 

1.413 

1.413 

.9984 

1.031 

*Note  here  the  asynptotic  nature  of  the  estimate  Eg  of  the  lowest 
achieveable  error. 
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Table  3 

Problem  A with  p = 2,  q = 1,  r = -1,  a = 1/100 
dhiform  mesh,  jIIUqUj,  = 0.28678 


m 

E 

^0 

6 

CO 

5 

25.215 

9.619 

.3270 

3.227(+3) 

10 

13.696 

4.809 

.4467 

4.640(+3) 

20 

7.296 

2.404 

.5966 

8.316(+3) 

40 

3.813 

1.202 

.7610 

1.055(+4) 

80 

1.959 

0.601 

.8947 

1.308(+4) 

Table  4 

Problem  A with  p = 2,  q = 1,  r = -1,  o = V 00 
/^proximately  optiml  mesh,  jIIUqHj,  = 0.286( 


m 

E 

^0 

0 

(0 

5 

9.994 

9.619 

.7679 

3.355 

10 

4,809 

4.809 

.9331 

1.666 

20 

2.409 

2.404 

.9828 

1.184 

40 

1.203 

1.202 

.9956 

1.049 

80 

0.601 

0.601 

.9989 

1.012 
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Table  5 


A:  Approxinately  optimal  mesh  for  m = 10 
: p = 0,  q = 1,  r = -1/4,  a = 1/100 
: p=2,  q=l,  r = - ./4,  a = 1/100 


j 

W 

Case  A1 

Case  A2 

0 

.0000 

.0000 

1 

. 00207 

.0127 

2 

. 00487 

.0340 

3 

.00877 

.0676 

4 

. 01443 

.1170 

5 

.02308 

.1862 

6 

.03732 

.2798 

7 

.06318 

.4025 

8 

.11781 

.5598 

9 

.26831 

.7569 

10 

1,00000 

1.0000 
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Table  6 


Problem  A with  p = 0,  q = 1,  r = -1/4,  a = 1/100 

Partitions  obtained  f^rom  the  uniform  mesh  by  successively  subdividing 
# the  first  interval  into  half 


m 

E 

Eo 

e 

0) 

10 

73.768 

11.306 

.2950 

2.34(+7) 

11 

58.853 

10.278 

.4705 

5.85(+6) 

12 

42.286 

9.421 

.6922 

1.45(+6) 

13 

27.357 

8.696 

.8457 

3.45(+5) 

14 

17,634 

8.075 

.9340 

7.49(+4) 

15 

13.580 

7.530 

.9447 

1.41(+4) 
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Table  7 


Problem  B with  a = 1,  p = 5/2 
iMform  mesh,  jllUgllj,  = .071070 


m 

E 

^0 

e 

Cl) 

5 

43.462 

32.317 

.9759 

1.126(+2) 

10 

22.080 

16.158 

.9939 

1.757(+2) 

20 

11.083 

8.079 

.9984 

7.568(+2) 

40 

5.547 

4.039 

.99924 

3.142(+3) 

80 

2.774 

2.019 

.99990 

1.281(+4) 

Table  8 

Problem  B with  a = 1,  p = 5/2 
/^proximately  optinel  mesh,  jllugllj.  = .071070 

m E Eg 

e 

CO 

5 

33.869 

32.317 

.9466 

1.577 

10 

16.519 

16.158 

.9694 

1.676 

20 

8.153 

8.079 

.9823 

1.755 

40 

4.049 

4.039 

.9894 

1.788 

80 

2.018 

2.019 

.9933 

2.437 
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Table  9 


Problem  B with  a = 5,  p = 9/10 
Uhifoim  mesh 


m 

E 

*^0 

e 

CO 

5 

49.477 

18.174 

.9059 

4.049(+3) 

10 

26.554 

9.087 

.9742 

1.229(+4) 

20 

13.530 

4.543 

.9934 

4.621(+4) 

40 

6.797 

2.271 

.9983 

1.808(+5) 

80 

3.403 

1.135 

.9995 

7.173(+5) 

Problem  B with  a = 5,  p “ 9/10 
Approxinately  optimal  mesh 

m E 

Table  10 

^0 

6 

CO 

5 

17.021 

18.174 

.7988 

2.617 

10 

9.181 

9.087 

.9217 

2.822 

20 

4.521 

4.543 

.9595 

2.324 

40 

2.254 

2.271 

.9820 

1.661 

80 

1.138 

1.135 

.9958 

1.614 
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Table  9 


Problem  B willi  a = 5,  p = 9/10 
Uniform  mesh 


m 

E 

=0 

e 

0) 

5 

49.477 

18.174 

.9059 

4.049(+3) 

10 

26.554 

9.087 

.9742 

1.229(+4) 

20 

13.530 

4.543 

.9934 

4.621(+4) 

40 

6.797 

2.271 

.9983 

1.808(+5) 

80 

3.403 

1.135 

.9995 

7.173(+5) 

Problem  B with  a = 5,  p = 9/10 
/^proxinately  optiiel  mesh 

m E 

Table  10 

Eq 

e 

0) 

5 

17.021 

18.174 

.7988 

2.617 

10 

9.181 

9.087 

.9217 

2.822 

20 

4.521 

4.543 

.9595 

2.324 

40 

2.254 

2.271 

.9820 

1.661 

80 

1.138 

1.135 

.9958 

1.614 
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•feble  U 


Problem  B:  /^proximately  optinal  mesh  for  m = 10 
Case  Bl:  a = 1,  p = 5/2 
Ccise  B2:  a = 5,  p = 9/10 


j 

Case  Bl 

Case  B2 

0 

.0000 

.0000 

1 

.0:37 

.4192 

2 

.1859 

.6918 

3 

.3001 

.7715 

4 

.5218 

.8255 

5 

.6872 

.8673 

6 

.7754 

.9016 

7 

.8442 

.9309 

8 

.9025 

.9565 

9 

.9538 

.9794 

10 

1.0000 

1.0000 
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